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Abstract 

We are concerned with a nonstandard phase field model of Cahn-Hilliard type. 
The model, which was introduced by Podio-Guidugli (Ric. Mat. 2006), describes 
two-species phase segregation and consists of a system of two highly nonlinearly 
coupled PDEs. It has been recently investigated by Colli, Gilardi, Podio-Guidugli, 
and Sprekels in a series of papers: see, in particular, SIAM J. Appl. Math. 2011 
and Boll. Unione Mat. Ital. 2012. In the latter contribution, the authors can treat 
the very general case in which the diffusivity coefficient of the parabolic PDE is 
allowed to depend nonlinearly on both variables. In the same framework, this pa- 
per investigates the asymptotic limit of the solutions to the initial-boundary value 
problems as the diff'usion coefficient a in the equation governing the evolution of 
the order parameter tends to zero. We prove that such a limit actually exists and 
solves the limit problem, which couples a nonlinear PDE of parabolic type with an 
ODE accounting for the phase dynamics. In the case of a constant diffusivity, we 
are able to show uniqueness and to improve the regularity of the solution. 
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Vanishing diffusion limit in nonstandard phase field systems 



1 Introduction 



In this paper, we consider the following system 



{l + 2g{p))^t^^ + ^^g'{p)^tp 
dtp-aAp + np) = pg'{p) 
p)V/i) ■ ulr = and 
/i(0) = po and p(0) = po, 



d„p\T = 



div(fi:(p, p) Vp) = 



(1.1) 
(1.2) 
(1.3) 
(1.4) 



in the unknown fields p and p, where the partial differential equations f ll.ip - fll.2p are 
meant to hold in a three-dimensional bounded domain Q, endowed with a smooth bound- 
ary r, and in some time interval (0, T). Relations fll.4p specify the initial conditions for p 
and p, while f ll.3p are nothing but homogeneous boundary conditions of Neumann type, 
involving precisely those boundary operators that match the elliptic differential operators 



This system has been recently addressed in the paper [6] : the existence of solutions 
has been proved, thus complementing and extending the results of the papers |3l IH [5] 
concerned with simpler or reduced versions of the problem. 

Here, we are interested to investigate the asymptotic behavior of the above initial- 
boundary value problem ( ll.ip - fll.4l) as the positive diffusion coefficient a appearing in 
([L2]) tends to 0. 

Let us briefiy explain the modelling background for fll.ip -f lT^p . Such a system comes 
from a generalization of the phase-field model of viscous Cahn-Hilliard type originally 
proposed in [Tl], and it aims to describe the phase segregation of two species (atoms 
and vacancies, say) on a lattice in presence of diffusion. The state variables are the 
order parameter p, interpreted as the volume density of one of the two species, and the 
chemical potential p. For physical reasons, p is required to be nonnegative, while the 
phase parameter p must of course take values in the domain of /'. 

We also recall the features of [3] and what has been generalized in [5], [6]. Firstly, the 
nonlinearity / considered in [3] is a double-well potential defined in (0, 1), whose derivative 
/' diverges at the endpoints p = and p = 1: e.g., for / = /i + /2 with /2 smooth, one 
can take 



with c a positive constant. In this paper, we let /i : M — [0, +oo] be a convex, proper 
and lower semicontinuous function so that its sub differential (and not the derivative) is 
a maximal monotone graph from R to M. Then, we rewrite equation (11. 2p as a differ- 
ential inclusion, in which the derivative of the convex part fi of / is replaced by the 
subdifferential (3 := dfi, i.e.. 



Note that /i need not be differentiable in its domain, so that its possibly nonsmooth and 
multivalued subdifferential /3 := dfi appears in (11.21) in place of /{. In general, /3 is only 
a graph, not necessarily a function, and it may include vertical (and horizontal) lines, as 
for example when 



in dnD-dig). 



flip) 



c (p log(p) + (1 - p) log(l - p)). 



(1.5) 



dtp-aAp + ^ + f^{p)= pg'ip) with ee/3(p). 



(1.6) 




(1.7) 
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and P = 5/[o,i] is specified by 

{< if p = 
= if < p < 1 . (1.8) 
> if p = 1 

Secondly, while in [3] (jf was simply taken as the identity map g{p) = p, in g is 

allowed be any nonnegative smooth function, defined (at least) in the domain where /i 
and its subdifferential live. The presence of such a function g allows for a more general 
behavior of (the related term in) the free energy, which reads 

^(p, Vp, p) = -| -pgip) + /(p) + || Vp|^ (1.9) 

Indeed, in particular g{p) is not obliged, as it was instead for g{p) = p, to take its minimum 
value at p = 0, be increasing and with maximum value at p = 1 (when -D(/i) = [0, 1]), 
but we may have many other instances like, e.g., a specular behavior of g around the 
extremal points of the domain of /. Here, we have to impose an additional restriction 
on g, which however looks reasonable from the modelling point of view: we postulate that 
g is a. (smooth) concave function, which in turn implies convexity with respect to p of 
the term —pg{p) in the free energy (11. 9p . However, let us recall that / may stand for 
a multi-well potential in which the nonconvex perturbations are incorporated into /2, so 
that in its entirety needs not be convex with respect to p. 

An important generalization that is considered in this paper concerns the diffusivity k. 
In [3], K was just assumed to be a constant function, but it can be a positive- valued, 
continuous, bounded, and nonlinear function of p (and this was the setting of [5]), or 
of p and p as it is postulated in |6]. For simplicity, we confine ourselves to study of 
the convergence properties of the solution under an assumption that guarantees uniform 
parabolicity, i.e., k, > > 0. We point out that [5] treats the situation of n depending 
only on p and possibly degenerating somewhere. 

Therefore, the system 



{1 + 2g{p)) dtp + fig'{p) dtp- div {^ip,p)Vfi)=0 (1.10) 

dtp-aAp + ^ + f^ip)=pg'{p) with ee/3(p), (1.11) 

(k(p, p)Vp) ■ i^lr = and d^p\r = (1-12) 

p(0)=po and p(0) = po, (1.13) 



turns out the initial and boundary value problem for a nonstandard and highly nonlinear 
phase field system in which however the role usually played by the temperature is here 
conducted by the chemical potential p. In the study of phase field systems, it has been 
always considered rather important to analyze the behavior of the problem as the coef- 
ficient a of the diffusion term in the phase parameter equation tends to 0. The limiting 
case a = corresponds indeed to a pointwise ordinary differential equation (or inclusion) 

5tP + e + /2(p) = /^^?'(p), ee/3(p), (1.14) 



in place of fll.lll) . and to an expression for the free energy f 1 1.9 1) in which the last quadratic 
term accounting for nonlocal interactions is removed. 
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In fact, especially for the choice f ll-Tp -f TTTS]) . the limiting problem can be formulated 
in terms of hysteresis operators: in particular, the so-called stop and play operators are 
involved; the interested reader can find some discussion and various results on this class 
of problems in[7l[8l|9l[IOl[IIl[T2l[l3]. 

By collecting a number of estimates independent of a for the solution {ficr,Pcr) to the 
problem f ll.lOp -f Tl.lSp . by weak and weak star compactness we prove that any limit in a 
suitable topology of a convergent subsequence of {(/^o-. Pa)} yields a solution to the limiting 
problem in which f ll.lip is replaced by f ll.l4p . Furthermore, under natural compatibility 
conditions on the nonlinearities and the initial data, we show boundedness for all the 
components of any solution to the limit problem. Finally, in the special case of a constant 
mobility k, in fll.lOp . we prove that the solution is unique and more regular than required. 

The paper is organized as follows. In the next section, we state precise assumptions 
along with our results. The basic a priori estimates independent of a are proved in Sec- 
tion |3] and they allow us to pass to the limit by compactness and monotonicity techniques. 
Finally, the last section is devoted to the study of the limit problem and our boundedness, 
uniqueness, and further regularity properties are proved. 



2 Assumptions and results 

The aim of this section is to introduce precise assumptions on the data for the mathe- 
matical problem under investigation, and establish our main result. We assume to be a 
bounded connected open set in with smooth boundary F (treating lower-dimensional 
cases would require only minor changes) and let T G (0, +oo) stand for a final time. We 
introduce the spaces 

V:=H\n), H:=L\n), W := {v e (Vl) : d^v = on T} (2.1) 

and endow them with their standard norms, for which we use a self-explanatory notation 
like II • II y. For powers of these spaces, norms are denoted by the same symbols. We remark 
that the embeddings W <Z V <Z H are compact, because VL is bounded and smooth. The 
symbol ( ■ , ■ ) denotes the duality product between V* , the dual space of V , and V itself. 
Moreover, for p G [1, +oo], we write || ■ ||p for the usual norm in LP(f2); as no confusion 
can arise, the symbol || • ||p is used for the norm in Lp{Q) as well, where Q := Q x (0,T). 

Now, we present the structural assumptions we make. It is useful to fix an upper 
bound for a, that is, 

0<CT<1. (2.2) 
Then, for the diffusivity coefficient k, we assume that 

K : (m, r) H- ^(m, r) is continuous from [0, +oo) x R to R, (2-3) 

the partial derivatives and d^n exist and are continuous, (2-4) 

K*, K* G (0,+oo), (2.5) 

< K{m,r) < K*, \drK{m,r)\ < k*, \d^K{m,r)\ < k* for m > and r G M, (2.6) 
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and for other nonlinearities we require that 

/ = /i + /2, /i [0,+oo], f,:R^R, (2.7) 

/i is convex, proper, l.s.c. and /2 is a function, (2.8) 

g G C2(M), g(r) > and g"{r) < for r e M, (2.9) 

/g, g, and g' are Lipschitz continuous. (2-10) 

It is convenient to introduce the notations 

k' := drK, K,":=d^K, (5:=dfi, and ir := (2-11) 

pm pm pm 

K{m,r) := / K(s,r)ds, Ki{m,r) := / K'(s,r)ds, K2{'m,r) := / K"(s,r)ds 
Jo Jo Jo 

for m > and r G M. (2.12) 

We denote by -D(/i) and -D(/3) the effective domains of /i and /3, respectively. Thanks to 
(12. 6p . it is clear that 

max{ I i^(m, r) I, 1/^1(771, r)|, \K2ifn,r)\} < tCm for every m > and r G M. (2.13) 

We also note that the structural assumptions of [5] are fulfilled if k only depends on m, 
and that, due to the presence of /3(p), a strong singularity in equation (11.111) is allowed. 
On the other hand, equation (II.IUI) is uniformly parabolic, since g is nonnegative and k 
is bounded away from zero. 

Remark 2.1. Let us recall that any convex, proper, l.s.c. function is bounded from below 
by an affine function (cf., e.g., [H Prop. 2.1, p. 51]), whence the assumption /i > looks 
reasonable, as one can suitably modify the smooth perturbation f2- Moreover, we point 
out that the sign conditions > and g" < are just needed on the set D{I3), for g can 
be extended outside of D{f3) accordingly. 

Concerning the initial data, we require that 

/io G V, /io > a.e. in Q, (2.14) 
PoeV, PoeDU\) a.e. infi, h{p^) e L\^) (2.15) 

and point out that the above assumptions regard the initial data for the limiting problem, 
i.e., the one with (11.141) in place of (II. lip . On the other hand, let us consider a family of 
initial data Pq^, poo- with 

/ioa e 1/ n L°^(fi), poa > a.e. in Q, (2.16) 
Poa G W, there is G H such that poo- G -D(/3), ^oa- G /3(poa-) a.e. in fi, (2.17) 

that approximate po, po in the sense that 

Poa Po and poa Po weakly in V, (2.18) 
||/i(pocr)||i is bounded independently of a. (2.19) 

For the reader's convenience, we show that such a family {poa, Poa} actually exists. Of 
course, if po ^ L°°{Q) we can take as poo- some truncation of po, e.g., poo- = minjpo, V^}- 
Concerning po„, one possible choice is letting poo- G W denote the solution to 

Poa - cr^poa + cr^oa = Po, wlth ^ (^{poa), a.e. in n. (2.20) 
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Indeed, the elliptic problem fl2.20p has a unique solution for all a > 0, since —A + /3 is a 
maximal monotone graph in H x H with effective domain 

{v eW : 3 r] e H such that v G D{P), t] G P{v) a.e. in Q}. 

Thus, poo- is nothing but the outcome of the application of the resolvent of —A + (3 to po 
(let us refer to [1] and [2] for basic definitions and properties of maximal monotone oper- 
ators). A formal test of the equality in (12.201) by ^o^t and the definition of subdifferential 
lead us to the estimate 

/ /i(po<x) + ^||Wlli< / /i(Po), (2.21) 
Jn Jn 

which ensures (12.171) and (12.191) . thanks to the nonnegativity of /i. A rigorous way of 
proving the existence of poo- and estimate (I2.2ip passes through the use of the Yosida 
approximation (see, e.g., [21 p. 28]) in place of (3. 

Now, we recall the result proved in [6] that allows us to specify a solution to the 
problem (I1.10I) - (I1.12I) . with a > 0, which fulfills the appropriate initial conditions. 

Proposition 2.2. Assume that both fl^ - (l2T^ and (I^T^ - d^TTD hold. Then, there 



exists at least one triplet {pa,Pa,^a) satisfying 

p^ G H^i'°°(0, T; H) n H\0, T; V) f] L°°(0, T; W), (2.22) 

GL-(0,r;/7), (2.23) 

p„ e H\0,T;H)nL°°{0,T;V)nL°^{Q), p„ > a.e. in Q, (2.24) 

div(K(/io-, Pcr)V/io-) G L^((5) and (^(/io-, Po-)V/i) ■ = a.e.onT,, (2.25) 

and solving the system of equations and conditions in the following strong form 

[l + 2g{p„))dtiia + fJ'a g'iPa) dtPa - div{K{fi^,p„)Vfia) = a.e.inQ, (2.26) 

dtPa - cr^Pa + + '^{Pa) = l^a9'{Pa) and ^„ E f3{p„) a.e.inQ, (2.27) 

Po-(O) = pocr and Pcr(O) = poo- a.e. in Q. (2.28) 



Let us point out that equation (I2.26P can be rewritten as 

dtU^ - div(K(p^,p^)V/i^) = Pag' (pa) dtpa, 

where = {I + 2g{pa))pa, a.e. in Q, (2.29) 



and the auxiliary variable has been added. Now, we take advantage of a variational 
formulation of (12.290 which also accounts for the boundary condition in (12.250 . that is, 

{dtUa{t),v)+ / (h;(p^,p^)Vp^)(t) ■ V?; = / Pa g\Pa) dtPaV 

Jn Jn 

for allv eV and a.a. t G (0, T). (2.30) 

The main result of this paper reads as follows. 
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Theorem 2.3. Assume that f l^ -l ET^ and (HH^ -I KT^ hold. For any a e (0,1] let 
ifJ'a, P(7,(,a) be the triplet defined by Proposition \2.1^ and let u„ := (1 + 2g{p„))ii„. Then, 
there exists a subsequence, still labelled by the parameter a, and a quadruplet (/i, p, u) 
such that 



weakly star zn L'^{0,T;H)r] L'^{0,T;V), (2.31) 

pa ^ p weakly star in H\0,T;H)nL°^{0,T;V), (2.32) 

-^^ weakly m L'^iQ), (2.33) 

u^^u weakly mW^^^/^{0,T;V*)r\L\0,T;W^'^/\n)) (2.34) 

as a \0. Moreover, any quadruplet (/i, p, ^, u) that is found as limit of converging subse- 
quences yields a solution to the following limit problem 

{dtu{t),v) + J^K{p,p)Vfi{t) ■ Vv = J^pg'{p)dtpv 

for allv eV and a.a. te (0, T), (2.35) 

u = {1 + 2g{p))fi a.e.inQ, (2.36) 

dtP + ^ + T^ip) = 1^ g'ip) and ^ e /3{p) a.e.inQ, (2.37) 

/i(0) = po and p(0) = po ci.e. in fl. (2.38) 



Remark 2.4. The nonnegativity property /i > a.e. in Q plainly follows from fl2.24p 
and dMH). 

Remark 2.5. One standard situation for the limit problem fl2.35l) - fl2.38l) is obtained for 
/3 = 9/[o,i] (cf. ([II7D-([IIHD). In this case becomes 

- 7r(p) + p g\p) - dtp e 9/[o,i] (p) a.e. in Q. (2.39) 

Then, if one introduces the generalized "freezing index" 

w{x,t) := / {-n{p) + iig\p)){x,s)ds, {x,t) e Q, 
Jo 

we thus have dtw — dtp G (9/[o,i](p), or equivalently, p = Sx[w^], where Sk is the stop 
hysteresis operator associated with the closed convex set K = [0, 1] (see, e.g., [IDl [HI [12] ) • 
Hence, we may rewrite fl2.39p as 

dtW = —ti^SkIw]) + pg^Sxlw]) a.e. in Q. 

In addition to the convergence result stated in Theorem 12.31 one can derive bound- 
edness for both the components p and ^ of any solution to the limit problem, provided 
that special additional requirements are satisfied, namely, by assuming that there exist 
real constants p*, p*, ^* such that 



p.,p*eD{f3), e*e/3(p*), re/3(p*), 

e* + 7r(p*)<0, r + vr(p*)>0, 

g'{p*)>0, g'{p*)<0. 



(2.40) 
(2.41) 
(2.42) 
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Theorem 2.6. In addition to the assumptions of Theorem \2.'^ suppose that ( I2.40p - fl2.42p 
and 

P* ^ Po ^ P* (i-^- iiT- ^ (2.43) 
hold. Then, the components p and ^ of any solution {fi,p,C,,u) to problem (12.351) - fl2.38p 
satisfy 

P*<P<P* and ^*<^<r a.e.inQ. (2.44) 

// moreover 

/io G L°°(fi) (2.45) 

and K = kq is constant, then the solution of Problem (12. 35p - (12.381) is unique and 

p G H\0, T- H) n L°°(0, T; V) n L'^ifd, T; W). (2.46) 

Remark 2.7. We observe that the above resuh is very general. Indeed, assumptions 
(I2.40p - (l2.42p are fulfilled with suitable constants for any graph /3 with bounded domain 
that generalizes the examples (II. 5p or (II. 7p . Of course, the decreasing function g' (cf. (12. 9p ) 
should not assume a definite sign on D{f3). 

Now, we list a number of tools and notations we owe to throughout the paper. We 
repeatedly use the elementary Young inequalities 

a6<7a^ + — 6^ and ah < -da^ + (1 - ■d)b^ 
47 

for every a, 6 > 0, 7 > 0, and ?9 G (0, 1) (2.47) 

as well as the Holder and Sobolev inequalities. The precise form of the latter we use is 
the following 

W^'P{n) C L«(fi) and \\v\\y < Cpjv\\w^,p^n) for every v G W^'P{n), 

provided that 1 < p < 3 and 1 < q < p* '■= (2.48) 

3 — p 

with a constant Cp^g in (I2.48P depending only on Q, p, and q, since Q C M.^. Moreover 

the embedding W^'P{Q) C L«(fi) is compact if l<q<p*. (2.49) 
The particular case p = 2 of (I2.48P becomes 

V C L\Q) and \\v\\q < C\\v\\v for every G and g G [1, 6] (2.50) 
where C depends only on Q. Moreover, the compactness inequality 

IK'llg < £^||Vi;||2 + Cg,e||t;||2 for every G g G [1,6), and e > (2.51) 

holds for some constant Cg^^ depending on Q, q, and e, only. We also recall the interpo- 
lation inequalities, which hold for any i) G [0, 1], 

ll^llr < ll^llj II^IIJ"'^ Vt; G LP{Q) n L^{Q), 

1 ?9 1-^9 

where p,q,r E [l,-|-oo] and - = — I . (2.52) 

r p q 

IK'||L'-i(0,T;L'-2(n)) < \\v\\LPi(^o,T;LP2{n)) h'\\]jii{0,T;Ln(n)) 

Vt; G LP'{0,T-LP^{Q))nL^'{0,T;L''^{Q)), 

1 1^ 1-1^ 

where , , rj G [1, +00] and — = 1 for z = 1,2. (2.53) 

ri Pi qi 
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We observe that fl^3^ implies \\v\\r < ^v\\p + (1 - ^)\\v\\g for every v G L^^Q) fl 
thanks to the Young inequahty, and a similar remark holds for fl2.53p . Thus, we have the 
continuous embeddings 

LP{Q)nL''{Q) CV{Q) and LP'{0,T; LP^{Q)) n L'''{0,T; L'^iQ)) C L"^(0,T; ^"^(fi)). 

We stress the important case of the space L°°(0, T; L'^{^)) nL^(0, T; L^{Q)), which occurs 
several times in the sequel and corresponds to pi = oo, p2 = 2, qi = 2, and q2 = 6. In 
particular, the choices 'd = 2/5 and = 1/7 yield the inequalities (for every v of the above 
space) and the continuous embeddings 

II^IIlio/3(q) < ll^llx^ll^^lly^ and X nY C L^^/^Q) (2.54) 
h'\\LVHo,T;L^^/Hn)) <h\\T\MT and X nY C L'^/^{0,T; L^^/'^{n)) (2.55) 
where X := L°°{0,T; L^{n)) and Y := L^{0,T; L%n)). 

Notice that we can take v G L°°(0, T; H) n ^^(0, T; V") in (l234D - fl235D . since V C L%n). 
Finally, we set 

Qt:=nx (0, t) for t G [0, T], (2.56) 

and, again throughout the paper, we use a small-case italic c for different constants, that 
may only depend on Q, the final time T, the shape of the nonlinearities / and g, and 
the properties of the data involved in the statements at hand; a notation like signals 
a constant that depends also on the parameter e. The reader should keep in mind that 
the meaning of c and q might change from line to line and even in the same chain of 
inequalities, whereas those constants we need to refer to are always denoted by capital 
letters, just hke C in fl230|) . 



3 The asymptotic analysis 

In this section, we prove Theorem 12. 3[ which ensures the existence of a solution to prob- 
lem fl2.35p - fl2.38p along with the convergence properties stated in fl2.3ip - fl2.34p . 

Then, for any a G (0, 1] we let (/io-, p^, C,a) denote the triplet defined by Proposition 12.21 
and set u„ := {l + 2g{pf^))fj,fj. The existence of (/ig-, p„, has been proved in [6]: we follow 
in parts the arguments developed there in order to recover useful estimates independent 
of a. Before that, let us remark that the property /io- > can be verified by simply 
multiplying equation 02.260 by — the negative part of Pa-, and integrate over Qf. In 
principle, in this computation one has to define k everywhere, e.g., by taking an even 
extension R with respect to the first variable. We observe that 

[(1 + 2g{p,{t))) dtp. + p. g{p.) dtp.] i-p-) = 1 dt{{l + 2g{p.m l/i^p) . 

Hence, by using Pq. > and owing to the boundary condition in fl2.25p . we have 

I [{l + 2gip.{t)))\p~{t)\^+ [ Rip.,p.)\Vp~\'' = fora.a. tG (0,T). 
^ Jn JQt 
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As both g and R are nonnegative, this imphes /i^ = 0, that is, /io- > a.e. in Q. 
First a priori estimate. We test (12.261) by /ig- and point out that 



[(1 + 2g{p^)) dtfi, + /i, g\p^) dtp.]fi^ = ^dt [(1 + 2g{p„)fil] . (3.1) 
Thus, by integrating over (0,t), where t G [0,T] is arbitrary, we obtain 

{l + 2g{p^{t)))M)\' + 2 I K{p,{s),p,{s))\Vp^^ = [ (l + 2(7(po<x))/iL. 
! JQt Jn 

We recall that g is nonnegative and Lipschitz continuous (cf. (I2.9p -( !2?T01) ). Moreover, 
Poa, Poa are both uniformly bounded in V by (12.181) . whence 

/ (1 + 2^(poa))^L < C {WpOaWl + WpOcrhWpOaWl) < C 

Jn 

owing to the Holder and Sobolev inequalities (see (12.501) ). Then, in view of g > and 
K > > 0, from (13.11) it follows that 

\\Pa\\L°°(0,T;H) + \\Pa\\L^{0,T;V) < C. (3.2) 

Second a priori estimate. We add pa to both sides of ( 12.270 and test by dtPa- On 
account of (I22D-(E8]) and (12111) . we obtain 

/ \dtp.\'+l\\pAt)rH+^\\'^pAt)rH+ 1 fiip.it)) 

= ^ / |Vpo.P+ / /(po<x) + ^ / (p^(t)-2/2(p.(t)))+ / p.^?V)5*p. 

^ Jo. ^ Jo. JQt 

for every t G [0,T]. Then, thanks to the Lipschitz continuity of and g, and owing to 
the bounds entailed by (I2.18p - (l2.19p . we find out that 

|5*P.P + ^llp.Wlli + f livp.(t)||i+ I WAt)) 



n 



<C + C I IPaWr + ^Z \dtPa\^ + c\\n4loo(^o^T;H)- 

Jn JQt 

On the other hand, by the chain rule and the Young inequality (I2.47p we have that 



c / IPaWP < c / |po<7 



|2 



\ [ \dtP.\' + c f\\pAs)\\lds. 

^ JQt JO 



n Jn 4 jQ 

Then, as /i is nonnegative, by accounting for ( 13.20 . with the help of the Gronwall lemma 
we infer that 

[ \dtp^\^ + \\pAt)\\l + a\\VpAt)\\l<c forall tG [0,T]. 

JQt 
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Thus, we conclude that 

\\Pu\\m{Q,T;H) + (^^^'^\\Pa\\L^{0,T;V) < C. (3.3) 



Third a priori estimate. We proceed formally and test (12.271) by —Ap„. Hence, 
integrating by parts and with respect to time, we deduce that 



^l|Vp.(t)||^ 



jQt JQt 

< \ I iVpoal'- / 7r'(p,)|Vp,p+ / (7'(p.)V/i, ■ Vp, + / /(p^KlVp^p, (3.4) 
^ Jn JQt JQt JQt 

where the equality ,^0- = (3{pa) has been used along with the smoothness of f3, according to 
our formal procedure. In fact, what is important is that the related term on the left-hand 
side is nonnegative, i.e., 

/ /3'(p.)|Vp.P>0. 
JQt 

Concerning the right-hand side of (13. 4p . we have that - f |Vpoo-|^ < c due to (12.181) . 

2 Jn 

and the estimate 



/ 7r'(p,)|Vp.P+ / ^?'(Pa)V/i.- Vp, < c / ||Vp.(s)||^ds + c||/i^||i2(o,r; 
JQt JQt JO 



V) 



owing to the boundedness of tt' and g' (see (I2.10I) - (I2.11I) ). About the last term, (12.91) and 
(12:241) imply 



JQt 



SO that the sign properties of g" and p^ become crucial to control this term. Then, in 
view of (13. 2p . from (13. 4p it follows that 

^ II Vpa(t) Hi + a / \Ap^\^ <c + c [ \\Vp^{s)\\lds for all te [0,T], 

2 JQt Jo 

and the Gronwall lemma and (13. 3p allow us to deduce that 

||P<7||L°°(0,T;y) + (^^^'^\\Pa\\L2(0,T;W) < C (3.5) 

Note that here we have used the regularity theory for elliptic equations, owing to the 
bound on (T||Apo-||| and to the homogeneous Neumann boundary condition satisfied by Po- 
(cf. (I2.22P ). Finally, an easy consequence of (13.31) and (13.51) comes out from a comparison 
of terms in (I2.27P , which yields 

||^a||L2(0,T;/f) < C. (3.6) 

Fourth a priori estimate. As Uo- = (1 + 2g{p(j))pai by (I2.10p we have that 

l^i^l < c(l + Ip^I) |/i<^|, 
IVm^I = \2g'{p,)p,Vpa + (1 + 2^(p,))Vp<,| < c | Vp,| + c (1 + \p,\) |V/i,|. 
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Now, taking fl3.2p into account, we see that |V/io-| is bounded in L^(0, T; L^(f2)), while 
is bounded in L^(0,T; L^{Q)) thanks to the Sobolev inequahty (12 .50 p . On the other 
hand, (jSlD provides a bound for | Vp^| in L~(0, T; ^^(fi)) and for \p„\ in L~(0, T; L^{n)). 
Hence, using Holder inequality, it is not difficult to check that the products |Vpcr| 
and |po-| |V/io-| are bounded in L^(0, T; L^/^(r2)), whereas \pa-\ \pa\ is even bounded in 
L^(0, T; Therefore, we conclude that 

ll'"'T||L2(o,r;iyi.3/2(f7)) < c. (3.7) 



Fifth a priori estimate. Let us recall that (13. 2 p and (I2.50p imply the boundedness 
of {pa} in the space L°°{0,T; L^{n)) n L^{0,T; L^{n)). Then, we can apply f l232|) with 
p = 2, g = 6, -i? = 1/2, r = 3 to see that 



IPCT 



< 



IP'C 



2||/i^(t)||6 for a.a. t G (0,T), 



whence squaring and integrating with respect to t lead to 

II Il4 II Il2 II Il2 

\\p(T\\L-iiO,T;L3(n)) ^ \\P'AL°°iO,T;L2{n)) II /^f^ II L2(0,T;L6(n)) ^ 



(3.^ 



Consider now (12.301) which turns out to be a variational formulation of (12.261) . As we 
want to prove that 

||(?t^ia||L4/3(o,T;y*) < c, (3.9) 

we use (12.301) and let v vary in L'^[0,T]V). By integrating with respect to time and 
invoking (12. 6p . the boundedness of g' and Holder's inequality, we obtain 



{dtu^it),v{t))dt 



< K*\\Vp^\\l^o,T;H)\\^v\\l^O,T;H) + C / || /i<^(t) || 3 1| (9tP^ (t) || 2 || I' (t) || 6 • 

Hence, in view of (13. 2p . by applying the Holder and Sobolev inequalities (see (I2.50p ) in 
the time integral, we infer that 

T 

{dtu,{t),v{t))dt 

< c||t'||L2(o_T;y) + c\\fia\\L-i{0,T;L'^n))\\dtPa\\L^{0,T;H)\\v\\L^{0,T;V)- 

Now, the continuous embedding L'^{0,T;V) C L'^{0,T]V), (13. 8p and (13. 3p allow us to 
conclude that 



{dtu^{t),vit))dt 



< c\\v\ 



L^{0,T;V), 



whence (13. 9p follows. 



Passage to the limit. By the above estimates, there are a quadruplet {fi, p,C,,u), with 
yLi > a.e. in Q, and a function k such that (I2.31I) - (I2.34I) are satisfied as long as 



K(/io-,Po-) — ^ k weakly star in L°°{Q) 



(3.10) 
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at least for a subsequence r = rj\0. By the weak convergence of time derivatives, the 
Cauchy conditions (12.281) hold for the limit pair (p, u). By (12.321) . (I2.34p . and the compact 
embedding (12.491) . we can apply well-known strong compactness results (see, e.g., [151 
Sect. 8, Cor. 4]) and infer that (possibly taking another subsequence) 

pa ^ P strongly in C°([0,T]; Lp{Q)) for p < 6 and a.e. in Q (3.11) 
Ua-^u strongly in L'^{0, T; Lp{Q)) for p < 3 and a.e. in Q. (3.12) 

The weak convergence (12.331) . together with (13. lip with p = 2, implies that C, G f3{p) a.e. 
in Q (see, e.g., [21 Prop. 2.5, p. 27]), due to the maximal monotonicity of the operator 
induced by /3 on L'^{Q). Now, we deal with the other nonlinear terms and the products. 
We first observe that (13.110 also entails that 

0(Pa) ^ <P{p) strongly in C0([0, T]; LP(fi)) for p < 6 and a.e. in Q (3.13) 

for = (yf, vr, 1/(1 + 2g), thanks to the Lipschitz continuity of such functions. This 
is sufficient to establish equation (12.371) . Indeed, by accounting for (I2.3ip . we see that 
the product Po-g{Pa) converges to pg{p) weakly (e.g.) in L'^iQ). On the other hand, 
(13.51) implies that crApa converges to zero strongly in L^{Q). Now, we prove equations 
(I2.35p - (l2.36p . which involve the whole triplet {p,p,u). The first step is showing strong 
convergence for p^ and relation (12.361) . By combining (I3.13P with (I3.12p . we see that 

u u 

Uo- = — T r-: a.e. in Q. (3.14) 

l + 2gipa) l + 2gip) ^ ^ ' 

This and (I2.3ip imply p = u/{l + 2g{p)) and (I2.36P is proved. Moreover, as {pa} is 
bounded in L^^/^{Q) by ([32]), the Sobolev embedding V C L^{n), and ([234]), we can also 
deduce a strong convergence. We summarize as follows: 

Pa ^ p strongly in L^iQ) for every p < 10/3 and a.e. in Q. (3.15) 

From this, we immediately infer that n{pa^Pa) converges to K,{p, p) a.e. in Q, just by 
continuity. Then, (I3.10p imphes k = hi{p, p) and 

t^{pa, Pa) ti{p, p) strongly in U'{Q) for every p < +oo. (3.16) 

Therefore, hi{pa, Pa)'^Pa converges to K{p,p)'Vp weakly in L^iQ) for every p < 2, thanks 
to (12.310 . and the choice p = 3/2 yields 

/ K(/i^, Pa)Vpa -Vv ^ I k{p, p)V p ■ Vv for every v G ^^(0, T; W^'^{n)). 
Jq Jq 

On the other hand, Pag'{Pa)dtPa converges to pg'{p)dtp weakly at least in L^{Q), as one 
can easily see by combining (12.321) . (13.130 . and (I3.15p . It follows that 

/ f^ag'{Pa)dtPaV / pg'{p)dtpv for cvcry G L~(Q). 
Jq Jq 

Moreover, (12.341) holds. Hence, we can conclude that 

{dtu(t),v{t)) dt + / n{p,p)Vp-Vv= / pg'{p)dtpv 
Jq Jq 

for every v G L3(0, T; W^'\n)) n L°^(Q). (3.17) 
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Now, we observe that dtu E L^/''^{0,T;V*) by (EMD and that k{p, p)Vv E L^{0,T;H) 
by flOTj) and the boundedness of k. Finally, pg'{p)dtp E L^/^{0,T; L'^/^{n)), since g' is 
bounded, dtp E L\0,T;H), and p E L\0,T;L%n)) as a consequence of flOTj) . V C 
L^(fi), and f l3.8p ). Therefore, we can improve fl3.17p by a density argument and see that 
the variational equation still holds for any v E L^(0,T; V). What we obtain is equivalent 
to f l2.35p . and the proof is complete. 

4 Properties of the limit problem 

In this section, we prove Theorem 12.61 In the whole section, it is understood that the 
assumptions of Theorem 12.61 are satisfied, and sometimes we do not remind the reader 
about that. As far as the first part of Theorem 12.61 is concerned, the true result regards 
ordinary variational inequalities and we present it in the form of a lemma. For convenience, 
we use the same notation p, etc., even though it is clear that everything is independent 
of x: the dot over the variable p denotes the (time) derivative, here. 

Lemma 4.1. Let fl2.40p - fl2.42p hold and p* < Po < p* ■ Then for every nonnegative 
function p E L^(0,T), the differential inclusion 

p{t) + f3{p{t)) + 71 {p{t))-p{t)g'{p{t)) 3 for a.a. t E {0,T) and p(0) = po (4.1) 

has a unique solution p E W^'^{0,T) such that 

P* < pit) < p* and < C(^) < r for a.a. t E (0, T), (4.2) 

where 

at) := -{p{t) + 7r{p{t))-p{t)g'{p{t))) E ^{p{t)). 

Moreover, there exists a constant C > such that if pi, p2 £ -^^(0,7) and pi, pi are two 
inputs and pi,p2 are the corresponding solutions of (14. ip . then for every t E [0,T] we 
have 

\Pl - P2\{t) + / \pl - p2\{T)dT 

Jo 

< C (ipS - Pol + ^ ((1 + Pi)IPi -P2\ + \pi- /i2|) (r) dr^ . (4.3) 

Proof. The existence of a unique solution can easily be proved, e.g., by the iterated Banach 
Contraction Principle, due to the monotonicity of /3 and to the Lipschitz continuity of 
the other nonlinearities. In (14. 2p . we only prove the upper inequalities since the proof of 
the lower ones is quite similar. It suffices to prove the desired inequalities for the solution 
(p, a of the cut-off problem 

p(t)+e(t) + vr*(p(t))-p(t)^*(p(t)) = 0, e(t) e /3(p(t)) fora.a. te(0,r), (4.4) 
p(0)=po, (4.5) 

where vr* and g* are defined by 

7r*(r) := 7r(min{r, p*}) and g*{r) := g'{mm{r, p*}) . 



Colli — Gilardi — Krejci — Sprekels 



15 



We test dm]) by (p - p*)+ and integrate. Recalling fCT?]) - (IT12]) and noting that 
and g*{p) = g\p*) where p > p*, we obtain 

lMt)-pT\' < -J\^-C){p-PT-J\c + ^*{P*)){P'PT 

+ [\n{p*)-n{p)){p-pT+ I\9\p){p-P*y 
Jo Jo 

< f (7r(p*) - n{p)) (p - p*y < c f I (p - p*)+p 



and the assertion is obtained by the Gronwall argument. The second inequality follows 
from the monotonicity of (5. Moreover, the lower bounds can be checked in a similar way. 
To prove (|13D, we set Wi{t) = fii{t)g'{pi{t)) - 7r{pi{t)), C.i{t) = Wi{t) - pi{t), i = 1,2. 
We have (.^i — ^2)(pi — P2) > almost everywhere. The function sign(,^i — ^2) (with 
sign(O) = 0) is bounded and measurable, and so is sign(pi — P2). We now claim that by 
testing the identity 

(6 - 6) + (Pi - P2) = Wi-W2 (4.6) 
by sign(^i — ^2), we infer that 

1^1 - 61 + 4:IPi - P2I < ki - W2I a.e. in(0,T). (4.7) 
at 

Indeed, this is obvious for all t such that sign(^i — 6)(t) = sign(pi — p2)(t) or such that 
^i(t) = i2{t). The remaining case is sign(^i — 6)(t) 7^ 0, sign(pi — p2)(^) = 0. For almost 
all t with this property, we have pi{t) = P2(i), ^|pi — P2|(^) =0, and (14. 7p follows. Using 
the Lipschitz continuity properties in (I2.10p and integrating ( 14. 7p over (0,t), we obtain 
for t G (0,T) 

^ \^i-^2\{s)ds+\pi-p2m <c(^\pI- pI\ + ((l + ^^)|p^_p2| + |pi-p2|)(r)rfr^ . 

On the other hand, ()4.6I) yields 

/ \pi - P2\{s) ds < / (|wi - W2I + 16 -6l)(s)c?s 
Jo Jo 

and (14.31) follows from the sum of the last two inequalities. □ 

Next, if (/i, p, 6 w) is a solution to problem (I2.35p -( l2.38p . it is clear that, for almost all 
X G fi, the functions p{x, ■) and p{x, •), and the constant po{x) satisfy the assumptions of 
Lemma 14. 1[ Thus, the first part of Theorem 12.61 concerning bounds (I2.44p is proved. We 
derive an interesting consequence. 

Corollary 4.2. Under the assumptions of Theorem \2.6[ let {p,p,$,,u) be a solution to 
problem (12. 35p - (12.381) satisfying the regularity conditions specified in Theorem \2.3[ Then 

peL°°{Q) and dtpeL°°{Q). (4.8) 



16 



Vanishing diffusion limit in nonstandard phase field systems 



Proof. We already know that both ^ and 7r(p) are bounded. Moreover, fig'{p) belongs 
to L°°{0,T;H) n L^{0,T;L^{n)) since /i does so and g'{p) is bounded. We see that 
also dtp belongs to such a space, just by comparison in (12.371) . It follows that dtp G 
L''/'^(0, T; L^^/'^(f2)) by fl2.55p . From this and assumption f l2.45p . we derive the bounded- 
ness of p. Indeed, we can reproduce the proof carried out in [B], Fifth a priori estimate], 
since that proof acts only on the equation for p and works provided that an estimate of 
dtp in L^/^(0,T; L^'^^^{^1)) is known. At this point, by comparing in (12.371) once more, we 
conclude that dtp is bounded as well. □ 

Remark 4.3. The analogous estimate 

P*< Pa< P* a.e. in Q (4.9) 
for the solution to problem (I2.26p - (l2.28p also holds provided that 

P* < Poa < P* a.e. in Q. (4.10) 

We prove one of the inequalities (14.90 . the other one being similar. We proceed as in the 
proof of Lemma [4.1^ testing (I2.27P by {p„ — p*)~^ and integrating. By accounting for the 
second inequality (I4.10p . we easily obtain 

I [ \{p.-pTit)\' + ^ [ \v{p^-pT\' 

^ Jn JQt 

+ [ {^.-C){p.-pT+ [ {C + n{p*)){p.-pT 

JQt JQt 
<[ (vr(p*)-7r(p.))(p.-p*)++ / p^g'{p^){p^-p*y. 



Now, we observe that all the terms on the left-hand side are nonnegative, the third one 
thanks to (I2.40p and the monotonicity of /3 (as before, the integrand vanishes whenever 
Per < p*), the last one due to (12.411) . Concerning the right-hand side, we show that the last 
integrand is nonpositive. Indeed, g' is decreasing (see (12.90 ). whence g\pa) < g'{p*) < if 
Pct > p*, and Per > 0. By taking all this into account and owing to the Lipschitz continuity 
of vr (cf. (12. lip ), we can apply the Gronwall lemma and conclude that {p„ — p*)^ = 0, i.e., 
p < p* a.e. in Q. 

Remark 4.4. A sufficient condition for (I4.10p to hold at least for small a is that poo- is 
given by (12.200 and the hypotheses of Theorem 12.61 are reinforced by also assuming that 

either infesspo>p* and supesspo<p* or ^* < < (4-11) 

The proof is rather simple and we show just one of the desired inequalities since the other 
one is quite similar. We test (I2.20p by (poo- — p*)^. We easily obtain 

/ \ipo.-pT\'+^ [ iv(po. -prr+^ / (w-r)(po.-pr 

Jn Jn Jn 

ipo-p*-^niPoa-pT- (4.12) 
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In the first case fl4.1ip . we set 5 := p* — sup ess po and take a* > such that a* |^*| < 5. 
Then, for a < a*, we have Po — p* — (t^* < —6 + a*\^*\ < a.e. in Q, so that the right- 
hand side of (14.121) is nonpositive. In the second case (14. lip , the same conclusion trivially 
holds. As the last two terms on the left-hand side are nonnegative (since (12.401) holds, 
/3 is monotone, and the third integrand vanishes whenever poo- < P*), we conclude that 
{Poa - p*)+ = 0, whence po^ < P*- 



Proof of the second part of Theorem 12.61 Assume thus that K,{fi, p) = kq and set 
for simplicity = I. The system now reads 

{dtu{t),v)+ I Vp(t)-Vi;= / fig'{p)dtpv 
Jn Jn 

for aAlv eV and a.a. t G (0, T), (4.13) 

u= {1 + 2g{p))fi a.e. in Q, (4.14) 

dtp + i + -n^p) = pig'{p) and ^ G /3(p) a.e. in (4-15) 

p(0) = Po and p(0) = po a.e. in f2. (4.16) 

Let (pi, Pi, ^i, Mj), = 1,2 be two solutions of (I4.13p -( l4.16p . We integrate (14.130 in time 
from to t and subtract the equation with index 2 from the one with index 1. We test 
the result hy v = (pi — P2)(i) and obtain, by virtue of Corollary 14. 2[ that 



{ui - U2)(pi - P2){t) + 



1 d 
2~dt 



V(pi - p2)(i 



r 



< c 



j {Ipi - P2\{t) (|pi - P2I + |pi - P2I + \dtpi - dtp2\) (r)rfr^ . (4.17) 



In addition, from Lemma [4.11 (see, in particular, (14. 3p ) and Holder's inequality it follows 
that 



n \Jo 



\dtpi - dtp2\{r) dr ] <c 



n \Jo 



t ^ 2 

- P2I + \Pi -P2|)(^) dr 



I \pi- P2?{s) <D I I {\pi- p2\^ +\pi- p2?){r)dT 
Jn Jo Jn 



(4.18) 
(4.19) 



for every t, s G [0,T], thanks to the boundedness for pi ensured by Corollary 14. 2[ Note 
that the constant D in (I4.19P is marked for later reference. 

Now, we observe that the inequalities 

(ni -n2)(pi -P2) > |pi-p2|^-2pi(5((pi)-5((p2))(pi-p2) > ^ |pi - P2r - c|pi - p2|^ 

hold a.e. in Q. Thus, by integrating (I4.17P from to s, s G (0, T), and ignoring a positive 
term on the left-hand side, we obtain 




\pi — p2\'^{t) dt < c 



Jn 




\Pi - P2p(^) dt + c 



Jn 




Jn 



Pi - P2p(i) dt 
1/2 



1/2 



X 



(4.20) 
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Hence, using Young's inequality and fl4.18p . we have that 




l/ii — /i2p(^) dt < c 



Jn 




\pi - P2?{t) dt 



Jn 



+c 




Jn \Jo 



x2 

ilP'i ~ f^2\ + \ Pi - P2\) {r)dr ] dt. 



(4.21) 



We now multiply (14.211) by 2D and add it to (I4.19p . Thus, we obtain an inequality of the 
form <l>(s) < c ^{t)dt, with 



$(s) 



\Pi - P2\ [s) 




l/ii - fJ'2\^{t) dt. 



Jn 



From the Gronwall argument, it is straightforward to deduce that $(s) = for all s, 
hence, /xi = fi2, pi = P2, which implies uniqueness. 

The bound for dtp can be established in the following way. Assume first that 
/io G W. We extend fi by /io and p by po for t < 0. Then, equation (I4.13P then can be 
written as 



{dtu{t),v) 



V/x(t) ■ Vv 



ip{t) V for all f G y and a. a. t G (0, T), 



(4.22) 



where '0 is defined by il){t) = (^fig' {p)dtp) (t) for t > and tjj{t) = — AyUo for t < 0. We 
observe that G L°°{—T,T]H) thanks to Corollary 14.21 and to our assumption on /iq. 
Next, we integrate (14.221) in time from [t—h) to t for any fixed t G (0, T) and a small h > 0, 
with the intention to let h tend to zero, and test the resulting equality by fi{t) — fi{t — h). 
We obtain 



{u{t) - u{t - h)) {i2{t) -^{t-h))+ - j - 



d 



V p{r) dr 



t-h 



n \Jt~h 



iPir)dr) {^(t) - fi{t - h)) 



ipir) dr 



W)-pi{t-h)\^^ 

<\jjKt)-pi{t-h)\^ + ch 



H 



(4.23) 



Now, we recall that (I4.14p holds, that g is nonnegative and Lipschitz continuous, and that 
/i and dtp are bounded by Corollary 14. 2 [ Hence, we easily derive that 

{u{t) - u{t - h)) {pit) - pit - h)) 

> |/i(t) - /i(t -h)\^- 2p.it) \gipit)) - gipit - h))\ |/i(t) - /i(t - h)\ 

> \fiit) - fiit - h)\'^ - ch\fiit) - fiit - h)\ > ^\fiit)-fiit-h)\'^-ch^. 

Therefore, by integrating (I4.23p from to T, forgetting the nonnegative term that in- 
volves V/i, and rearranging, we obtain 




\fxit) - fxit ~ h)f dt < ch"^ + c 



Jn 



V/io dr 



<cK'. 
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As > is arbitrarily small, this implies that dtfi G L'^{Q). At this point, we are allowed 
to use the Leibniz rule for the time derivative dtu; then, from f l4.13p - fl4.14p we infer that 
the equation 

(1 + 2gip))^t^^ + fig'{p)dtp - = (4.24) 

holds at least in the sense of distributions. By comparison, we deduce that Afi G L'^{Q), 
whence G L^(0,T; W). Using the identity 

- J^dtfiAfi = ~JjVfi\' a.e. in(0,T), 

we see that V/i G L°°(0,T; ^^(fi)). Thus, the regularity (EJSD is estabhshed if /iq G W. 

Let now /iq G V fl L°°{Q) be arbitrary, and consider a sequence C W bounded 

in L°°{Q) and converging to /iq in as — oo. Let {^k, Pk,^k,Uk) be the corresponding 
solutions to fl4.13p -f l4.16p . Then, we can use equation (14.241) written with the index k and 
test it by dt^k- We obtain 

/ ^^'^'^^^^'^Ijtf l^'"'^'^^)!'^ / \Mt)\\^t^^km (4.25) 

with an obvious choice of tpk G L?'{Q) bounded in this space (even better) independently 
of k. By time integration, it is straightforward to obtain a bound for ||9t/ifc||L2(Q) and for 
II V/Xfc||L°o(o,T;_ff) independent of k. Then, by weak star compactness we infer that 

/ifc /i weakly star in H^{Q, T; H) n L°°(0, T; V) 

at least for a subsequence, which implies (see, e.g., [151 Cor. 4, p. 85]) strong convergence 
in C°([0, T]; if). In particular, /i(0) = /iQ. On the other hand, {fik, Pk,^k,Uk) satisfies the 
estimates stated in Lemma Wl] and the boundedness properties for /i^ and dtPk given by 
Corollary 14.21 which are uniform with respect to k. This yields weak or weak star limits 
p and C,- Moreover, strong convergence in L^{Q) for {pk} and {dtPk} is ensured via a 
Cauchy sequence argument based on (14. 3p . integration over Q, and Gronwall's lemma. 
Hence, {pk}, {Pk}, {dtPk} converge strongly in Lp{Q) for every p G [l,oo). At this 
point, it is not difficult to verify that {p, p, ^, u), with the corresponding u, actually solves 
problem fl2.35l) - fl2.38l) and thus coincides with the unique solution (/i,p, ^,m). Therefore, 
the proof is complete. 
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